Signal detection by means of phase coherence induced through phase resetting 
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Detection and location of moving prey utilizing electrosense or mechanosense is a strategy com- 
monly followed by animals which cannot rely on visual sense or hearing. In this article we consider 
the possibility to detect the source of a localized stimulus that travels along a chain of detectors at 
constant speed. The detectors are autonomous oscillators whose frequencies have a given natural 
spread. The detection mechanism is based on phase coherence which is built up by phase resetting 
induced by the passing stimulus. 

PACS numbers: 87.10.+e,87.19.Bb 



The ability to detect, locate, and capture prey is vital 
for survival. Many animals accomplish these tasks us- 
ing visual or acoustic information. However, species that 
have developed in an environment where these senses are 
obscured, have to rely on alternative mechanisms. For 
example, the paddlefish (Polyodon spathula), found in 
the river basins of the Midwestern United States and in 
the Yangtze River in China, makes use of a passive elec- 
trosensory system jl]. Another example is the weakly 
electric fish that combines active and passive electrosense 
with a mechanosensory lateral line system In these 
animals, receptors transform stimuli into electric signals 
which excite the terminals of primary afferent neurons. 
These afferents are well known to exhibit periodic spike 
patterns ||. 

In the last decades a lot of research has been devoted 
to the details of information processing on the neural 
level, i.e., the dynamics of single neurons or neural net- 
works. However, at the behavioral level still many open 
problems exist. Since the performance and the analysis 
of experiments usually involve an enormous effort, effi- 
cient and tractable models are indispensable, both for 
planning and interpretation. 

Here we present an idealized, however, analytically 
tractable model, proposing a mechanism for the detec- 
tion of a localized stimulus. This stimulus is passing an 
array of receptors, which we model as phase oscillators. 
To measure the degree of coherence between the oscilla- 
tors we choose the well known synchronization index [Q . 
First we examine the influence of a random initial dis- 
tribution of the oscillator phases on the synchronization 
index and introduce a threshold value to distinguish a 
stimulus from a "false alarm" . Then we investigate the 
influence of our model parameters for the detection of a 
moving stimulus. 

We consider a linear chain of N uncoupled phase 
rotors which are characterized by the set of variables 
ipi, . . . , tpjf- The rotors are aligned at equal distance 
along an axis of length L, i.e. the position of rotor k 
is x k = (N - k)/(N ~ 1)L, k = 1,...,N. Each rotor 



has its own natural frequency fi& which, in the absence 
of a stimulus, determines the simple linear growth of the 
phase, i.e. ipk(t) = V'fc(O) + &kt- We assume the fre- 
quencies to be independently and identically distributed 
according to a Gaussian with mean S7o and standard de- 
viation A£7 = t/Qq. 

An appropriate quantity to measure the degree of 
phase coherence among these rotors is the complex vari- 
able 
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This global order parameter contains both the informa- 
tion about the instantaneous collective phase <!>(£) and 
the instantaneous degree of phase coherence measured by 
the modulus R(t) at time t. Its square can be expressed 
in several ways: 
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This quantity is termed synchronization index since it is 
widely used in the description of synchronization pro- 
cesses Q]. From Eqs. (||) and (||) it is obvious that 
< R(t) < 1 with R(t) = 1 indicating perfect coher- 
ence. 

We initialize the array by randomly selecting a phase 
for each of the rotors according to the uniform distribu- 
tion on [0, 2%). Thus, the quantity Rq = R(t = 0) is a 
random variable. Its density contains important informa- 
tion because even in the absence of any signal the array 
of rotors will generate nonvanishing values of i?o- These 
have to be discriminated from values of R(t) which signif- 
icantly indicate coherence induced by the passing stim- 
ulus. Figure [l] shows numerically estimated densities, 
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FIG. 2: Sketch of the setup. The stimulus moves at constant 
speed v relative to the oscillator chain. Each time it passes 
an oscillator the phase is reset to ip* . 



FIG. 1: The density of the random variable 7?o for different 
N = 10, 100, 1000, 10000 (top down). The dots indicate the 
result of numerical simulations whereas the full line identifies 
the Rayleigh distributions ([?]) . Vertical lines mark the thresh- 
old set by the value 2/\fN, which corresponds to a 2% level 
of false alarm. 



where N = 10, 100, 1000, 10000 equidistributed phases 
were used to compute a single realization of the random 
variable R . An analytic expression for these distribu- 
tions can be derived by applying the central limit theo- 
rem (Lindenberg-Levy theorem) to the following pair of 
random variables 

1 N 1 N 

X ^^E cos [^(°)]< Y N = (5) 
fc=i fc=i 

which yields for large N the limiting density [H 

p(X N = x,Y N = y)~- exp [~N{x 2 + y 2 )] . (6) 

Changing to polar coordinates (Rq,$o) and integrating 
over the angle <J> immediately leads to the Rayleigh dis- 
tribution (Fig. [l]) 



Pn (Rq) = 2NR exjp(-NR 2 ). 
Mean and variance of this distribution read 
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a(R th ) = J 2NR exp(-NR 2 ) dR = exp(-NR 2 h ) 
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can be used to define a threshold value R th by demanding 
that a(i? th ), which is the probability for false alarm, be 



less than some fixed small number. For instance, i? th = 
2/VjV corresponds to a(R th ) < 2% which means values 
larger than 2/vN occur by random configuration with a 
probability of less than 2%. In what follows we will use 
R tb = 2/ y/N to discriminate stimuli against the random 
configuration background. 

A standard model describing phase resetting by an ex- 
ternal stimulus of strength / is given by the following 
phase dynamics H : 



ip = Cl + I cos ip . 



(10) 



This dynamics can be illustrated as the overdamped mo- 
tion in a tilted corrugated potential landscape. If I < 
no troughs (minima) and barriers (maxima) exist and 
the phase continues cycling forward (f2 > 0) at varying 
speed. For / > 51 two fixed points emerge, which corre- 
spond to a minimum at n — arccos(51/7) and a maximum 
at 7r + arccos(f2/I). For constant / the phase settles in 
the minimum (mod27r) regardless of the initial position, 
which means the phase eventually is reset to the corre- 
sponding value ip* — n — arccos(57/7). The situation is 
harder to analyze with a time varying stimulus I(t); the 
net effect will depend on many details of the stimulus, 
e.g., the time scale of variation, the height of the signal 
peak, etc. 

Our detection setup would require to consider N such 
phase equations each with its own time varying stimulus 
Ik(xo ~ %k ~ vt), where xq and v are the initial position 
(at time t = 0) and the constant velocity of the travel- 
ing stimulus, respectively. Irrespective of the details, the 
equation of motion will be too complicated to be solved 
analytically in closed form. If, however, the peak value of 
the stimulus is sufficiently high and the duration is short, 
we can simplify the resetting mechanism: The passing 
stimulus resets the phase tpk(t) to some global value ip* 
the very moment it is at position Xk, i.e., the reset is in- 
stantaneous. After this reset the phase again increases 
linearly with its natural frequency The situation is 
sketched in Fig. 0. The history of phase ipk can thus be 
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written as 



detector chain, namely, 



, M J MV + Ukt t<t k 

(for all k = 1, . . . , N), where t k is the time when the 
stimulus passes the oscillator k. Substituting this into 
Eq. (Q) we find the following value of the synchronization 
index: 

R 2 (t) = 4 + 172 i S kk(t) + S kN (t) + S NN (t)} (12) 
Jy iV z 

in the time interval t k <t < t k +\, where we denote 

k k 
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These expressions depend on the initial phases ipi (0) and 
the natural frequencies fij. We consider both quanti- 
ties to be random parameters of the model. To charac- 
terize the net effect of observing many realizations, i.e., 
to evaluate the mean performance of many individuals, 
we average the synchronization index over both the ini- 
tial phases (equidistributed) and the natural frequencies 
(Gaussian). The first average over the phases yields (for 
t k <t < t k+ i) 

k k 

(R 2 (t)) = ^ + J* £ £ cos " **) " - ^ ■ 
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Note that the value of tp* is irrelevant for this expression. 
The second average over the natural frequencies results 
in 



x exp 
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In the following we relate time to the position of the 
stimulus x(t), 



t = 



x a - x(t) 



(18) 



We can then derive an expression that reflects how the 
twice averaged global synchronization index varies as a 
function of the position of the stimulus over the linear 
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for x k +i < x < x k . Assuming the oscillators to be dis- 
tributed along the linear chain in an equidistant man- 
ner, i.e., x k = (N — k)Ax for k = 1, ...,7V with 
Ax = L/(N -1), we find 

fc(x) k(x) 
i=l j>i 

f v 2 * 2 A 

x exp<^ —g(x,i,j) } , (20) 



where 



Ax 



(21) 



and k(x) = rnin{int[./V — -£z],N\. The parameter n 
turns out to be related to the ratio of the travel time 
between two neighboring oscillators AT = Ax/v and the 
mean rotation period To = 2tt/Qq, i.e., n — 2tt^-. It is 
useful to write K — 2irm+5, where m 6 N and 8 E [0, 27r). 

Equations ( p0| ) and (21) present the central result of 
our model. The detection regions in the x-6 plane, 
i.e. where ((i? 2 (x/ £ ))) is larger than the threshold value 
i? 2 h = A/N, is shown in Fig. | for TV = 100 (top) and 
N = 10 (bottom). It can be seen that detection works 
only as long as detuning, quantified by 8, and frequency 
spread, coded by 77, are not too large. Moreover, we find 
that the detection region shrinks in the 8 direction, but 
enlarges in the x direction with increasing N, i.e., detec- 
tion already works when the stimulus has passed only a 
small number of oscillators. For small 77 and small 8 we 
can consider the following limiting cases: First let us deal 
with the case of zero frequency spread, i.e., 77 = 0. The 
double sum over cosines can be performed yielding the 
following expression, 



N-k + 
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(22) 



which we exemplify for N = 100 in Fig. ^. Depending 
on the detuning parameter S, constructive or destruc- 
tive effects show up. Introducing the frequency spread 
77 > erodes both the constructive and destructive ef- 
fects. Note that the cycle number m matters if r\ > 
whereas it is irrelevant for the case r\ = 0. In Fig. |5| we 
exemplify how the detection curve for TV = 100, m = 1 
and a detuning value of 8 — 0.017T is pushed below the 
detection threshold by an increasing frequency spread 77. 
These results indicate that the detection mechanism is 
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FIG. 3: Regions in which prey is detected in the x-S plane 
defined by the demand that ((R 2 (x))) > R^ = 4/N. Top 
panel for N = 100, r\ — (black) and r\ = 0.01 (checkered 
plotted on top). The bottom panel compares the region for 
N = 10 and ry = 0, 0.01, 0.02, 0.03 (light to dark gray plotted 
on top of each other) with N = 100 and r\ = 0.01 (checkered). 
Each time the stimulus passes an detector, ((R(x))} changes 
discontinously and decays for r) > 0. 
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FIG. 4: Variation of the global synchronization index as a 
function of the stimulus position x in the case of vanishing 
frequency spread ij = for N = 100. Depending on the 
detuning parameter S constructive or destructive effects of the 
array or rotors can be observed: 5/ir = 0, 0.01, 0.02, 0.03, 0.04 
shown as solid, dotted, dashed, long-dashed, dot-dashed lines, 
respectively. The solid horizontal line marks R? h — 4/N 
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FIG. 5: Spatial variation of the global synchronization index 
as a function of the stimulus position x in the case of varying 
frequency spread r\ for slight detuning, i.e., 8/ir — 0.01, and 
for m = 1 and N = 100. An increasing spread erodes the 
detection mechanism: 77 = 0, 0.005, 0.01, 0.015 shown as solid, 
dotted, dashed, long-dashed, respectively. The solid horizon- 
tal line marks R^ = 4/N 



rather sensitive with respect to the width of the frequency 
distribution for a large number ./V of oscillators. However, 
we would like to point out, that the biological relevance 
is not eradicated by this finding, since evolutionary opti- 
mization offers an explanation how the confined param- 
eter range might have been realized. 

In conclusion, we have presented a simplified but an- 
alytically tractable model for signal detection, which 
works by creating significant coherence in a chain of 
phase oscillators. This coherence is induced by a strongly 
localized stimulus that travels at constant speed and re- 
sets phases instantaneously. The ability to detect a stim- 
ulus rapidly is balanced by the sensitivity to variations 
in the oscillator frequencies or deviations from the opti- 
mal velocity. The variations in the frequencies, however, 
guarantee a fast desynchronization after the stimulus has 
passed. 

Although our approach concentrates on seemingly 
crude assumptions, it catches the main features of prey 
detection. Future experimental studies have to reveal in 
which direction this model has to be extended to account 
for given biological applications. 
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